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ABSTRACT 



Context. Formaldehyde is an organic molecule that is abundant in the interstellar medium. High deuterium fractionation is a common 
feature in low-mass star-forming regions. Observing several isotopologues of molecules is an excellent tool for understanding the 
formation paths of the molecules. 

Aims. We seek an understanding of how the various deuterated isotopologues of formaldehyde are formed in the dense regions of low- 
mass star formation. More specifically, we adress the question of how the very high deuteration levels (several orders of magnitude 
above the cosmic D/H ratio) can occur using H 2 CO data of the nearby p Oph A molecular cloud. 

Methods. From mapping observations of H2CO, HDCO, and D2CO, we have determined how the degree of deuterium fractionation 
changes over the central 3' x 3' region of p Oph A. The multi-transition data of the various H 2 CO isotopologues, as well as from 
other molecules (e.g., CH3OH and N 2 D + ) present in the observed bands, were analysed using both the standard type rotation diagram 
analysis and, in selected cases, a more elaborate method of solving the radiative transfer for optically thick emission. In addition to 
molecular column densities, the analysis also estimates the kinetic temperature and H2 density. 

Results. Toward the SMI core in p Oph A, the H 2 CO deuterium fractionation is very high. In fact, the observed D 2 CO/HDCO ratio is 
1.34 ± 0.19, while the HDCO/H 2 CO ratio is 0.107 ± 0.015. This is the first time, to our knowledge, that the D 2 CO/HDCO abundance 
ratio is observed to be greater than 1. The kinetic temperature is in the range 20-30 K in the cores of p Oph A, and the H 2 density 
is (6 - 10) x 10 5 cm 4 . We estimate that the total H 2 column density toward the deuterium peak is (1 - 4) x 10 23 cnT 2 . As depleted 
gas-phase chemistry is not adequate, we suggest that grain chemistry, possibly due to abstraction and exchange reactions along the 
reaction chain H 2 CO — > HDCO — > D 2 CO, is at work to produce the very high deuterium levels observed. 

Key words, astrochemistry - ISM: abundances - ISM: clouds - ISM: individual objects: p Oph A - ISM: molecules 



1. Introduction 

The study of deuterated molecules in the interstellar medium 
(ISM) has been intensified over the past decade ever since 
it was discovered that singly and multiply deuterated species 
(like the D-containing versions of H 2 CO, CH 3 OH, NH 3 ) oc- 
curred at abundances that were orders of magnitude higher than 
woul d be expected f rom the local ISM D/H ratio of about 1 .5 x 
10~ 5 dLinskvl 120031) . The regions that show these elevated D- 
abundances are mainly associated w ith low-mass pro t ostars . For 
insta nce, around IRAS16 293-2422, ICeccarelri et ail (119981) and 
later lLoinard et al.l (|2000) detected lines from doubly-deuterated 
formaldehyde (D?C O) which prev iously had only been s een in 
the Orion KL ridge dTurnerlll990l) . Later, lLoinard et at] d2002l) 
extended the study to a larger set of protostars arriving at 
D2CO/H2CO abundance ratios as high as 0.05-0. 4. Likewise, 
high degree s of de uterium fra ctionation in CH3O H (Parise et al 



2002 12004 120061) a nd NH 3 dRoueff et al.l l2000t lLoinard et al 
200UlLis etaTff 2002) were subsequently discovered 



These studies are all related to single telescope pointings. 
To our knowledge, only one effort to delineate the D2CO/H2CO 
abundance variat ion within a single sou rce (IRAS 16293-2422) 
has been made dCeccarelli et all 1 120011) . The D 2 CO/H 2 CO ra- 



* Based on observations with the Atacama Pathfinder Experiment 
(APEX) telescope. APEX is a collaboration between the Max-Planck- 
Institut fur Radioastronomie, the European Southern Observatory, and 
the Onsala Space Observatory 



tio peaks some 20" away from the protostar. Further out, the 
degree of deuteration seems to be lower in the quiscent gas 
dCeccarelli et alJl2002l). In the case of deuterated ammonia, the 
mapping study by Roueff et"al] (2005) revealed that the deu- 
terium peak is typically offset from the positions of the embed- 
ded protostars. These authors argued that the observed scenario 
could be explained by the formation of deuterium-enriched ices 
during the cold pre-collapse phase. At a later stage, a newly 
formed protostar evaporates the ices. 

In their extensi ve study of deuterated H2CO and CH3OH, 
iParise et all 12006) concluded that formation of CH 3 OH on 
grain surfaces was a likely explanation for the high degrees of 
deuterium fractionation seen in the various isotopologues. For 
H2CO, the situation was less clear and a formation path in- 
volving gas-phase reactions could not be ruled out. In fact, for 
the warmer Orion Bar PDR region. [Parise et"al 1 (120091) very re- 
cently advocated that gas-phase chemistry is entirely responsi- 
ble for the deuterium fractionation seen for some singly deuter- 
ated species (including HDCO). It should be noted that singly- 
deutera ted species h ave been detected in dark and translucent 
clouds dTurnerJl2001l) where the deuteration occurs only in the 
dense parts. 

In interstellar clouds, deuterium is mostly molecular in the 
form of HD. Deuterium can be transferred from this main molec- 
ular reservoir to other molecules via exothermic reactions with 
the molecular ions H3 , CHJ , and C2H2 . These basic reactions 
are followed by efficient ion-neutral reactions and, in some 
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cases, by reactions on grain surfaces. The exothermicity of the 
reactions in volving these three molecular ions are in the range 
230-550 K (iGerlich et al.l20^lAsvanvlTaT]l2004tlHerbst et all 
119871) . As a result, deuterium fractionation is efficient in cold 
environments. It is also well established that a high degree of 
depletion on the grains of CO, O, and other heavy species, 
which would otherwise destroy efficiently H^ (and its deuter- 
ated analogues), is another important condition for deuterium 
fractionations taking place. Correlations between CO depletion 
and fractionation have been obse rved towards prestellar cores by, 
e.g., iBacmann et al.l (120031) and ICrapsi et al.1 (l2005[). and in th e 
envelope of Class protostars by Emprechtinger et al. (2009). 
Because of these two conditions (low temperature and high CO 
depletion), deuterium fractionation is particularly efficient dur- 
ing the early stages of star formation. The surface reactions 
on cold grains important for deuterium fractionation of H2CO 
have recently been investigated through laboratory experiments 
dHidaka et al.ll2009h . 

The aforementioned low-mass protostar IRAS 16293-2422 
is located in the eastern part of the p Ophiuchi cloud complex. 
To the west in the s ame complex, mo re than 1 degree away, lies 
the p Oph A clo ud (Loren et al. 1990) at a distance of about 120 
pc dLombardi et al.l2008l:lLoinard et alj2008USnow et alj2008l) . 
This cloud core is well-studied by i nfrared, submillimeter, 
and millimeter continuum observations dWard-T hompso n et al] 
1989; lAndre et all Il993t iMotte et all 119981). It hosts a well- 
collimated molecular o utflow ( Andre et alj 1 1990h and its driv- 
ing source VLA 1623 dAndre et alJ 1993[). Moreov e r, it w as in 
the direction of this cloud core that lLarsson et alj |2007), us- 
ing the Odin satellite, detected the 119 GHz line from molecular 
oxygen as well as the ammonia ground state line at 572 GHz 
dLiseau et al.ll2003[). More recently , while searching for the 234 
GHz 18 line lLiseau et al.l (2010) detected several lines due to 
D2CO toward the millimeter continuum peaks. T his study also 
includ es C ls O(3-2) mapping observations. Earlier. lLoinard et all 
d2002l) reported a high D 2 CO/H 2 CO ratio toward the VLA 1623 
source. The existence of D2CO in several positions of this 2' x 3' 
cloud core formed the incentive of the present study as an ex- 
cellent source to delineate the distribution of H2CO, HDCO, and 
D 2 CO as we know the distribut ion of the dust dMotte et al.l 19981) 
and the gas in terms of C ls O dLiseau et al.ll2.Q10l) . Here we re- 
port on mapping observations of the p Oph A cloud core in sev- 
eral frequency settings that cover most of the low-energy H2CO, 
HDCO, and D2CO l ines in the 1 .3 milli meter band using the 12 
m APEX telescope (Gii sten et al.ll2006l) . In addition, lines from 
many other species were observed simultaneously (eg. CH3OH, 
SO, and N2D + ). Here the CH3OH results are of importance be- 
cause CH3OH is directly involved in the H2CO chemistry and 
the N2D + results are, of course, of interest for the deuteration. 
Although the sulphur chemistry is not of immediate interest here 
we chose to include the SO and SO2 observational and analysis 
results since they are of importance as a complementary tool for 
determining the physical conditions. This paper is organized as 
follows. In Sect. (f2]i we describe the observations and then, in 
Sect. ([3), we present the observational results. In Sect. (01 we 
obtain the physical conditions. Before concluding, in Sect. ([6j, 
we discuss our results in Sect. ©. 



2. Observations 

The APEX telescope at Chajnantor (Chile) was used to map the 
3' x 3' (with a spacing of 30") area centred on the coordinates 
a(J2000) = 16 h 26 m 27.2 s and £(J2000) = -24° 23' 34" which 
is very close to the SM1N position in the p Oph A cloud core 



as designated bv lAndre et al.l d 19931) in their submillimeter con- 
tinuum maps. We used the single-sideband tuned APEX-1 re- 
ceive r which is part of th e Swedish heterodyne facility instru- 
ment dVassilev et al.l200"8T) . It has a sideband rejection ratio more 
than 10 dB. The image band is 12 GHz above or below the ob- 
serving frequency depending on whether the tuning is optimized 
for operation in the upper or lower sidebands, respectively. For 
two frequency settings, the strong CO(2-l) and C 18 0(2-l) lines 
entered via the image band and from the strength we could es- 
timate the sideband rejection ratio to be about 15 dB in both 
cases. As backen d we used two 1 GHz modules of the FFTS 
dKlein et al.l2006l) . Each FFTS 1 GHz module has 8192 effective 
channels and the modules can be placed independently within 
the IF band width of 4-8 GHz. At 230 GHz this channel spacing 
corresponds to 0.16 kms~' and the telescope beamsize (HPBW) 
is 27". 



The telescope control software APECS dMuders et al.ll2006l) 
was used to control the raster mapping. All observations were 
performed in position-switching mode using a reference position 
offset by 300" east and 200" north of the map center. The point- 
ing of the telescope was maintained and checked regularly by 
means of small CO(2-l) cross maps of the relatively nearby car- 
bon stars IRAS 15 194-5 115 and RAFGL1922. The determined 
pointing offsets were generally consistent from day to day and 
we believe we have an absolute pointing uncertainty better than 
5". To optimize the antenna focussing we used continuum ob- 
servations on Jupiter and Saturn. 

The observations took place in two blocks and one additional 
day during 2009: April 24 - May 1, May 21, and July 4 - July 9. 
The column of precipitable water vapour was typically around 
0.7 mm (varied between 0.3-2.9 mm). Typical system tempera- 
tures for the frequencies in question (218-252 GHz) were 200- 
220 K. The telescope ma in beam efficiency is j] m b = 0.73 at 345 
GHz (iGiisten et al.ll2006h . Hence, using the antenna surface ac- 
curacy of 18 pm we estimate that (using the Ruze formula) that 
the APEX main beam efficiency around 230 GHz is just slightly 
higher than at 345 GHz, about 0.75 which we adopt when con- 
verting the observed intensities from the T* A scale to the T m \, 
intensity scale. The heterodyne calibration procedure at APEX 
is more elaborate than the standard chopper wheel calibration 
scheme and involves three measurements by apart from the nor- 
mal sky observation it also measures the receiver temperature by 
observing a hot and cold loads. Moreover, the atm ospheric con- 
tribution is based on the model dPardo et al.ll2001l) that has been 
adapted to the atmospheric characteristics at the Chajnantor site. 
The absolute calibration uncertainty is estimated to be 10% in 
the 1 mm band. 

In Table [1] we summarize the targeted formaldehyde lines. 
The additional lines are summarized in Table [2] In the tables 
we include the line parameters: the transition frequency (typi- 
cal uncertainty is 0.05 MHz or better), energy of lower level, 
and Einstein A-coefficient. These have been com piled from 
the Cologne Database of Molecular Spectroscopy (iMiiller et al.l 
I2001ll2005l) . We here also indicate the symmetry due to the nu- 
clear spin direction of the H (or D) nuclei, which for H2CO, 
Hj 3 CO and D2CO can be ortho (parallel spins) or para (anti- 
parallel spins). The statistical weight ratio of the symmetries is 
also noted if applicable. In the case of CH3OH, the internal rota- 
tion of the methyl group results in the symmetry species A and 
E. Note that radiative or non-reactive collisional transitions are 
forbidden between levels of different symmetry. 
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Table 1. Observed H 2 CO, H' 3 CO, HDCO, and D 2 CO lines 



Frequency 


Transition 


Ei 


A 


a 




Symmetry 


(MHz) 




(K) 


(s" 


i , 








H 2 CO (o/p 


= 3) 








218222.19 


3o,3 ~ 2 _2 


10.6 


2.8 x 


10" 


-4 


para 


218475.63 


32,2 — 2 2 ,1 


57.6 


1.6 x 


io- 


-4 


para 


218760.07 


32,1 - 22,0 


57.6 


1.6 x 


10 - 




para 


225697.78 


3i,2 - 2ii 


22.6 


2.8 x 


10- 




ortho 




H' 3 CO (o/p 


= 3) 








219908.52 


3i,2 - 2i,i 


22.4 
HDCO 


2.6 x 


io- 


-4 


ortho 


227668.45 


1 i.i _ Oo.o 


0.0 


1.7 x 


io- 


-6 




246924.60 


4l.4 - 3 U 


25.8 


4.0 x 


io- 


-4 






D 2 CO (o/p 


= 2) 








221191.66 


41.4-3,3 


21.3 


2.9 x 


io- 


-4 


para 


231410.23 


4o,4 - 3o,3 


16.8 


2.8 x 


io- 


-4 


ortho 


233650.44 


42.3 — 3 2 ,2 


38.4 


2.7 x 


io- 


-4 


ortho 


234293.36" 


4 3 . 2 - 3 3 ,i 


65.4 


1.6 x 


io- 


-4 


para 


234331.06° 


4 3 ,i - 3 3 ,o 


65.4 


1.6 x 


io- 


-4 


para 


245532.75 


4i. 3 - 3,.2 


23.1 


3.9 x 


io- 


-4 


para 


" only observed at (0, -30") 



Table 2. Additional lines 



Frequency 


Transition 


Ei 


A 






Symmetry 


(MHz) 




(K) 


(s- 


') 








CH 3 OH (A/E = 


1) 








218440.05 


42,2 — 3 1.2 


35.0 


4.7 x 


10- 


-5 


E 


241700.22 


5o,5 _ 4(),4 


36.3 


6.0 x 


10- 


-5 


E 


241767.22 


5l,5-4l, 4 


28.8 


5.8 x 


10- 


-5 


E 


241791.43 


5o,5 _ 4(),4 


23.2 


6.0 x 


10- 


-5 


A 


241879.07 


5l,4-4,.3 


44.3 


6.0 x 


10- 


-5 


E 


241904.15 


52,4 _ 42,3 


49.1 


5.1 x 


10- 


-5 


E 


241904.65 


52,3 — 4 2 .2 


45.5 
SO 


5.0 x 


10- 


-5 


E 


219949.44 


5 6 - 4 5 


24.4 


1.3 x 


10- 


-4 




251825.77° 


65 - 5 4 


38.6 
34 SO 


2.0 x 


10- 


-4 




246663.47 


65 - 54 


38.1 
SO, 


1.8 x 


10- 


-4 




226300.03 


14 3 .n - 14 2 .i2 


108.1 


1.1 X 


10- 


-4 




241615.80 


5 2 .4 - 4 1.3 


12.0 


8.5 x 


10- 


-5 




245563.42 


10 3 ,7 - 102,8 


60.9 

34 so 2 


1.2 x 


10- 


-4 




219355.01 


Hi, 11 - 10o,io 


49.5 


1.1 X 


10- 


-4 




241985.45 


83,5 - 82,6 


42.8 


1.1 X 


10- 


-4 




246686.12 


4 3 ,i - 4 2 , 2 


18.7 
N 2 D + 


8.4 x 


10- 


-5 




231321.83 


3-2 


11.1 


7.1 x 


10- 


-4 




" only observed at (0, -30") 



3. Results 

In this section we will first display the formaldehyde mapping re- 
sults of the p Oph A cloud core, both as spectra, integrated inten- 
sity maps or velocity position diagrams where appropriate. After 
that, the mapping results for the other molecules will presented. 
Comparisons will be made with the ex isting APEX C 18 0(3-2) 
data at 329 GHz of iLiseau et all d2010h and the IRAM 30 m 
continuum map at 1.3 mm of Mo tte et a The angular 

resolution of the two data sets is similar; the C 18 data have 
an HPBW of 19" while the 1.3 mm continuum map has an an- 
gular resolution of 15". Both these maps are shown in Fig. Q] 
where the 1.3 mm continuum map is shown in contours on top 
of a grey-scale image of the C 18 0(3-2) integrated intensity. The 




Aa [arcsec] 

Fig. 1. Colour image: The C 18 Q(3-2) integ r ated in tensity map 
of p Oph A cloud core from ILiseau et al J (120101) . The inten- 
sity scale is s hown to the right. Contours: The 1.3 mm contin- 
uum data bv iMotte etail (Il998l) . The contours start at 0.15 Jy 
with subsequent contours at every increment of 0. 15 Jy. The flux 
densities are given in a 15" be a m. W e show the positions of 
the mm peaks from IMotte et ail dl998l) and the C 18 peaks of 
ILiseau et al . (2010). The map offsets are given relative the posi- 
tion or(J2000) = 16 h 26 m 27.2 s and 5(J2000) = -24° 23' 34". 



3' x 3' region shown in Fig. [TJ corresponds to the area that has 
been mapped here in formaldehyde. 

3. 1 . Formaldehyde results 

All lines listed in Table[TJwere detected except the 1 i i -Oo,o tran- 
sition of HDCO at 227 GHz. This transition has a much lower 
spontaneous rate coefficient than the other transitions. Two of 
the D2CO lines were not mapped and were only observed in the 
(0", -30") position. The typical rms in a map spectrum is in the 
range 0.06-0.1 K. 

In Fig.|2]we plot the H 2 CO(3 , 3 - 2 0j2 ) and D 2 CO(4 ,4 - 3 0;3 ) 
spectra. They consist of 49 spectra separated by 30" on a 7x7 
grid. The H 2 CO lines posess a relatively complicated structure 
as compared to the D 2 CO lines. The latter line shows only one 
velocity component around 3.8kms _1 with full width at half 
maximum (FWHM) of 0.7kms _1 and it is peaking strongly at 
the (0", -30") position which, within the positional uncertain- 
ties, coincides with the SMI or P3 position (cf. Fig . [TJ . We will 
hereafter refer to this position as the D-peak. The D-peak is also 
clearly seen in the integrated intensity maps of the three different 
H 2 CO isotopologues (Fig.0. 

iMaret et alJ d2004l) present a H 2 CO(3 , 3 - 2 , 2 ) IRAM 30 m 
spectrum toward the position of VLA1623. The nearest posi- 
tion to VLA1623 in our map is at (0",-60"). The shape of 
our spectrum at th is position is similar to the one presented by 
IMaret et all d2004t) and our T m \, peak temperature of 3.8 K is 
slightly higher than theirs T*/rj m b ss 2.2/0.62 « 3.5 K despite 
the higher resolution of 11" in the IRAM 30 m spectrum. This 
would be expected for a relatively extended source which shows 
little variation on the scale of 1 1" - 27". Also, from the distribu- 
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p Oph A H 2 C0 and D 2 C0 
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Fig. 2. H2CO(3o,3 - 2o,2) and D2CO(4o,4 - 3o,3) map spectra toward the p Oph A cloud. The (0", 0") position is the same as in Fig.Q] 
The vertical scale in each spectrum represents T m b scale in K and the horizontal scale is velocity vlsr with respect to local standard 
of rest in km s _1 as shown in the upper rightmost panel. 



tion of the H2CO emission we do not see any strong component 
that can be attributed to VLA16230 

The velocity as obtained from the NqH+ fl-O) observations 
bv lDi Francesco et all (120041 12009|) (see also lAndre et~aT]|2007l) 
toward SMI is 3.7 km s _I with an FWHM of 0.6 km s -1 , i.e. very 
similar to our values determined from the D2CO lines at the D- 
peak. The peak intensity of the D2CO line is about half that of 
H2CO. The size of D2CO 4q 4 - 3o,3 emission region (see Fig. [3}, 
is found, by fitting a 2-dimensional gaussian, to be 33" X 61". 
This size corresponds to a deconvolved source size of about S = 
19" x 55" for a beam size of B — 27". The filling factor when 
pointing the 27" beam toward the center of the source, is then 
S 2 /(S 2 + B 2 ) = (19" x 55")/(33" x 61") ~ 0.5. This filling 
factor is merely an upper limit because any unresolved small- 
scale dumpiness may decrease the filling factor further. 

In order to investigate the H2CO line distribution over the 
p Oph A cloud core we plot four velocity position diagrams 
(Fig. |H between AS = ±90" at Aa = 0". Here we also see 
that HDCO is peaking at the D-peak at A6 = -30". Looking 
at the HDCO emission, there appears to be a velocity gradient 
over the D-peak source of about 0.2kms _1 when going from 
A6 - 0" to Ad = -60". It is also close to the D-peak where the 
C 18 0(3-2) emission has its maximum. The secondary C 18 0(3-2) 



Based on the radial intensi t y dist ributions of the sub-millimeter 
continuum, Jayawardhana et al. (2001) suggested that the infall zone 
around VLA1623 is surrounded by a constant-density region which 
is the dominating contribution to the radial intensity profiles at scales 
> 10". 



peak (PI in Fig. [TJ is at Ad = +30" and with a lower velocity of 
w 3.1 kms~' as compared to the velocity of +3.7kms~ 1 at the 
D-peak position and it is also here where the low-energy H2CO 
line intensities reach their maximum value at vlsr = 3.3 km s~'. 
We denote this peak PI from now on. There is also weak emis- 
sion emanating from both HDCO and D 2 CO at the A6 = +30" 
position. 

Interestingly, there is a third H2CO peak at AS = 0" and 
vlsr ~ +2.9 km s -1 which has no obvious counterpart in the 
C ls O velocity position map. This peak makes the H2CO(3o,3 - 
2o,2) spectrum at (0",0") to look doubly peaked (Fig. 0. 
However, these are two different cloud components and the dip 
is not an effect of self-absorption. This is evident from the veloc- 
ity position diagrams in Fig. [4] but is also clear when comparing 
the H' 3 CO(3i, 2 - 2 U ) and H 2 CO(3 u - 2 U ) spectra. In Fig. El 
the central three spectra from the ortho lines H2CO(3i 2 - 2i i) 
and H^CO(3i, 2 - 2 U ) are displayed. Toward the (0",6") posi- 
tion we see the doubly peaked line profile also for this H2CO 
line, however there is no hint that the H} 3 CO line is peaking at 
the velocity of the dip for H2CO. If anything, the emission from 
the rarer species seems to follow that of the main species. We 
will return to the low-velocity component seen in the (0",0") 
position below when we present the CH3OH results. 

3.2. CH3OH results 

In Fig. [6] we show a map of the integrated intensity for the 
CH3OH 5i 5 - 4j 4 and 42 2 - 3i 2 E-lines. The latter line comes 
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Aa [arcsec] 



Fig. 3. Maps of the H 2 CO(3 , 3 - 2 0j2 ), HDCO(4i ;4 - 3 U ), and 
D2CO(4o,4 - 3o,3) integrated line intensity over the range 3.0 - 
4.5kms _1 . 



from levels of higher energy (Table Q]). The distribution of 
CH3OH is quite different from that of H 2 CO (see Fig. 0. It 
should be noted that the CH 3 OH 4 2 ,2 - 3 1,2 line and the H 2 CO 
lines at 218 GHz were observed simultaneously so the differ- 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




Aa [arcsec] 



Fig. 6. Integrated intensity maps of the CH3OH 5 1,5 - 4; 4 and 
42,2 _ 3 1_2 E-lines. First contour is at 0.3 Kkm s and the incre- 
ment is 0.6 Kkm s _1 . 



ent distributions cannot be a result of large pointing offsets. The 
CH3OH emission has its maximum at (-30", 0) and in Fig . [7] all 
observed 5-4 lines in this position are shown. Here the peak 
emission velocity is lower, around vlsr ~ 2.8 km s _1 . The cloud 
component at (-30", 0") extends into adjacent positions. The 
low-velocity feature seen for H2CO (Fig. HJ is very likely as- 
sociated with the CH30H-peak. The higher energy H2CO lines 
32,2/1 - 2 2 ,i/o peak near this position. In addition, there is little 
CH3OH at the D-peak position (Fig.©. 

3.3. Results from other molecules 

Two lines from SO(56 - 4$) and 34 SO(65 - 54) were covered 
during the observations. The map spectra are shown in Fig. [9] 
There is a single strong peak in the 34 SO line intensity toward 
the position (-60", +60") which we will call the S-peak. Here 
the 34 SO and SO2 lines are very narrow, only about 0.4 km s _1 . 
This is about twice the width one would expect from thermal 
broadening only (at 7\ in = 20 K). Further west, at (-90", +30"), 
the 34 SO line is broader, about 1.1 km s -1 . The profile of SO 
shows line wings here which could be related to outflow ac- 
tivity. However, it is at the very border of our map and a clear 
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p Oph A 



5 o.5- 4 o.t A 



5 0.5-* ,4 E 



( -30.0, +0.0) 



5? jyq 4, 3/, E 



Jo* 



Table 3. Properties of the H 2 CO, CH 3 OH, SO, and S0 2 emission 
in selected cloud positions from gaussian fits 



241.7 241.75 241.8 241.85 

Rest frequency [GHz] 



Fig. 7. CH 3 OH 5-4 lines around 241 GHz toward the (-30", 0") 
position. Note that the 5 2i 4 - 4 2j 3 and 5 2i 3 - 4 2 2 E-lines are 
blended. 



delineation into red and blue wings are impossible. It should 
be pointed out that the S-peak has no structural counterpart in 
the C ls O(3-2) map nor the 1.3 mm continuum map (cf. Fig.Q]). 
Toward the D-peak, the SO lines exhibit a similar emission ve- 
locity (vlsr ~ 3.7kms _1 ) and line width (Av w 0.83 km s~') as 
the other molecules. The covered S0 2 and 34 S0 2 lines (Tabled 
show a distribution similar to that of SO and 34 SO. At the Pl- 
position the SO and S0 2 emission velocity of about » 2.9 km s _1 
is lower than the velocity of m 3.3kms _1 for the formaldehyde 
isotopologues and CH3OH. 

The deuterated version of N 2 H + (ob served by 
iDi Francesco etal.ll2004l l2009t lAndre et alJ l2007h . N 2 D + (3-2) 
at 231 GHz was also covered during the formaldehyde obser- 
vations. The integrated intensity map of N 2 D + (3-2) is shown 
in Fig. [TO] The peak intensity of the N 2 D + emission is clearly 
associated with the D-peak, with a secondary weaker source 
to the south-east. This secondary peak is not coincident with 
SM2 (see Fig. [TJ but is located in the direction where the dust 
emission extends. Moreov er, it shows up in the N 2 H + map of 
IDi Francesco et all d2004l) and is designated N6 by them. The 
emission velocity of N 2 D + (3-2) at the D-peak is 3.7 km s" 1 and 
the FWHM is 0.85 km s~' in agreement with what is found for 
the other species. 



3.4. Result summary 

After presenting the results of the distribution of the different 
molecules toward the p Oph A cloud above we will now make 
a summary of the results which will be used in the analysis be- 
low. For H 2 CO and CH3OH we have identified three positions 
of interest which seem to represent distinctive sources; the D- 
peak at (0", -30"), the PI peak at (0", +30"), and the CH 3 OH 
peak at (-30", 0"). In addition to these three peaks, the sulphur 
containing molecules are localized to (-60", +60") (the S-peak). 
The emission velocities and FWHMs as obtained for the differ- 
ent species in these sources are summarized in Table [3] (from 
gaussian fits of the lowest energy lines in those cases where we 
have mapped multiple lines). Integrated line intensities for the 
relevant positions are tabulated in Table [U Table [5] and Table [6] 
In these tables the Icr uncertainties of the integrated intensities 
due to noise have been entered. 



Offset 


Molecule 


Vlsr 


Av (FWHM) 






( km s" 1 ) 


(kms" 1 ) 


(0", -30") 


H 2 CO 


3.8(0.1) 


0.95(0.05) 




HDCO 


3.6(0.1) 


0.70(0.05) 




D 2 CO 


3.8(0.1) 


0.68(0.02) 




CH3OH 


3.7(0.1) 


0.73(0.05) 




SO 


3.7(0.1) 


0.83(0.05) 




S0 2 


3.7(0.1) 


0.67(0.05) 


(0",+30") 


H 2 CO 


3.4(0.1) 


1.20(0.05) 




HDCO 


3.2(0.1) 


0.85(0.12) 




D 2 CO 


3.2(0.1) 


0.60(0.15) 




CH 3 OH 


3.4(0.1) 


1.18(0.05) 




SO 


3.0(0.1) 


1.06(0.05) 




so 2 


2.8(0.1) 


0.77(0.05) 


(-30", 0")° 


H 2 CO 


2.8(0.2) 


0.95(0.10) 






3.6(0.2) 


0.87(0.10) 




CH3OH 


2.6(0.2) 


0.62(0.10) 






3.2(0.2) 


1.27(0.10) 


(-60", +60") 


SO 


3.0(0.1) 


0.78 (0.05) 




34 SO 


2.9(0.1) 


0.43(0.02) 




so 2 


2.9(0.1) 


0.42(0.02) 




34 S0 2 


3.0(0.1) 


0.35(0.02) 



Fitted by two components, full intensity used in models 



Table 4. H 2 CO, H^CO, HDCO, and D 2 CO integrated line in- 
tensities in three positions of the p Oph A cloud 



Molecule 




Line 


Frequency /,„(, 


= J 


T mb dv 


A/ mb • 








(MHz) 


(Kkrns" 1 ) 


(Kkms 1 ) 








D-peak (0", -30") 






H 2 CO 


3(1.3 


- 2o,2 


218222 




6.07 


0.05 




3 2 .2 


-2a,! 


218476 




0.83 


0.05 




3 2 ,i 


- 2 2 ,o 


218760 




0.80 


0.05 




3|.2 


-2i,i 


225698 




6.88 


0.06 


H> 3 CO 


3 1,2 


-2i,i 


219908 




0.26 


0.06 


HDCO 


4l,4 


- 3l,3 


246924 




0.93 


0.08 


D 2 CO 


4m 


-3,, 3 


221192 




0.87 


0.08 




4(1.4 


— 3o,3 


231410 




2.21 


0.05 




4 2 ,3 


— 32 2 


233650 




0.45 


0.005 




4 3 ,2 


-3 3 j 


234293 




0.032 


0.005 




43.1 


- 3 3 ,o 


234331 




0.028 


0.005 




4,,3 


-3i, 2 


245533 
PI (0",+30") 




0.57 


0.05 


H 2 CO 


3(1.3 


— 2o,2 


218222 




8.69 


0.05 




3 2>2 


-22.1 


218476 




1.28 


0.06 




32.! 


- 22,0 


218760 




1.19 


0.06 




3 1,2 


-2i,i 


225698 




10.96 


0.07 


H' 3 CO 


3,. 2 


-2i,i 


219908 




0.26 


0.07 


HDCO 


1m 


-Oo,o 


227668 






0.06 




4l,4 


" 3l,3 


246924 




0.44 


0.05 


D 2 CO 


4, .4 


-3 13 


221192 




0.13 


0.05 




4 ,4 


- 3q,3 


231410 




0.19 


0.05 






CH 3 OH-peak(-30",0") 






H 2 CO 


3o,3 


- 2o,2 


218222 




8.56 


0.05 




3 2 , 2 


-22,1 


218476 




1.45 


0.06 




32.! 


- 22,0 


218760 




1.41 


0.06 




3.,2 


-2i,i 


225698 




10.13 


0.07 


H' 3 CO 


3.,2 


-2i,i 


219908 




0.27 


0.06 


HDCO 


4l,4 


-3i, 3 


246924 




0.12 


0.05 


D 2 CO 


4, .4 


- 3l,3 


221192 




0.12 


0.05 




4 ,4 


_ 3o,3 


231410 




0.15 


0.05 



lcr error due to noise only 
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N„D + 3-2 




Table 6. Integrated line intensities for other molecules in four 
positions of the p Oph A cloud 







Fig. 10. Integrated intensity map of the N2D + 3-2 line toward 
p Oph A. First contour is at 0.2Kkms _1 and the increment is 
O^Kkms- 1 . 

Table 5. CH3OH integrated line intensities in three positions of 
the p Oph A cloud 



Line 


Frequency 


Imb - J T m b dv 


A/ mb " 




(MHz) 


(Kkms- 1 ) 


(Kkms- 1 ) 







D-peak(0",-30") 






4 2 ,2- 


3i, 2 El 


218440 


0.21 


0.05 


5o,5 _ 


4o,4El 


241700 


0.13 


0.05 


5i,5- 


4i, 4 E2 


241767 


0.72 


0.05 


5o,5 - 


4o,4A 


241791 


0.81 


0.05 


5l,4 _ 


4i, 3 El 


241879 


0.09 


0.05 


52,3/4 


- 42,2/3 El 


241904 


0.07 


0.05 






PI (0", +30") 






4 2 ,2 - 


3i, 2 El 


218440 


0.61 


0.05 


5o,5 _ 


4o, 4 El 


241700 


0.59 


0.05 


5l,5 - 


4i, 4 E2 


241767 


2.39 


0.05 


5o,5 - 


4o,4A 


241791 


2.88 


0.05 


5 1.4 _ 


4i, 3 El 


241879 


0.20 


0.05 


52,3/4 


- 4 2 ,2/ 3 El 


241904 


0.31 


0.05 






CH 3 OH-peak (-30", 0") 




4 2 ,2 - 


3i, 2 El 


218440 


1.43 


0.05 


5o,5 - 


4o, 4 El 


241700 


1.49 


0.05 


5i,5- 


4i, 4 E2 


241767 


4.53 


0.05 


5o,5 ~~ 


4o, 4 A 


241791 


5.25 


0.05 


5l,4 - 


4i, 3 El 


241879 


0.96 


0.05 


52,3/4 


- 4 2 ,2/3El 


241904 


1.43 


0.05 



lcr error due to noise only 



4. Analysis 

For several of the detected species (Table |4| we have mul- 
tiple transitions that have significantly different energy lev- 
els. In such cases the so-c alled rotation diagram analysis (eg. 
iGoldsmith & Langerffl999h can be employed to determine the 
rotation temperature, r rot , and the molecular column density, 
N mo \. If all lines are optically thin, the rotation diagram method 
often is adequate to analyse multi-transition data. However, here 
it is likely that we have a mixture of optically thin and thick 
lines and therefore we adopt the modified approach to the rota- 
tion diagram method described by Numme lin et al.l (|2000). The 



Transition 


Frequency I mb = J T mb dv 
(MHz) (Kkrns" 1 ) 


A/ mb a 
(Kkms- 1 ) 




D-neak <0" -3f)"1 






SO(5 6 - 4 5 ) 


219949 


7.22 


0.06 


SO(6 5 - 5 4 ) 


251826 


3.37 


0.04 


34 SO(6 5 - 5 4 ) 


946663 


0.19 


0.03 


S0 2 (5 2 , 4 - 4 13 ) 


241616 


0.69 


0.02 


SO 2 (10 3 ,7 - io 2 , 8 ) 


745563 


0.16 


0.04 


N 2 D + (3 - 2) 


231322 


3.38 


0.06 




PI fO" +30") 






SO(5 6 - 4 S ) 


219949 


13.1 


0.06 


34 SO(6 5 - 5 4 ) 


746663 


0.55 


0.04 


S0 2 (5 2 , 4 - 4 U ) 


241616 


1.70 


0.05 


SO2OO37 - 10 2 , 8 ) 


245563 


0.57 


0.04 


N 2 D + (3 - 2) 


231322 


0.37 


0.08 




CH 3 OH-peak(-30",0") 






SO(5 6 - 4 5 ) 


219949 


8.25 


0.06 


34 SO(6 5 - 5 4 ) 


246663 


0.08 


0.04 


S0 2 (5 2 ,4 - 4 U ) 


241616 


0.23 


0.05 


SO 2 (10 37 - 10 2 , 8 ) 


245563 


0.05 


0.04 


N 2 D + (3 - 2) 


231322 


0.17 


0.07 




S-peak (-60", +60") 






SO(5 6 - 4 5 ) 


219949 


14.5 


0.06 


34 SO(6 5 - 5 4 ) 


246663 


1.31 


0.04 


S02(5 2 ,4 -4,, 3 ) 


241616 


3.76 


0.05 


SO 2 (10 3 , 7 - 10 2 , 8 ) 


245563 


1.14 


0.04 


S0 2 (14 3 ,„ - 14 2 , l2 ) 


226300 


0.27 


0.03 


34 S0 2 (4 3 ,, - 4 2 , 2 ) 


246686 


0.24 


0.04 


34 S0 2 (8 3 , 5 - 82,6) 


241985 


0.18 


0.04 


34 S0 2 (ll U i - 10 ,,o 


) 219355 


0.46 


0.04 


NoD + (3 - 2) 


231322 




0.07 



l<r error due to noise only 



modification involves the inclusion of the peak optical depth and 
thus, in addition to T rot and N mo \, also the beam filling factor, rj b f , 
can be determined by minimization of a^- 2 -value. For a gaussian 
source distribution, rj b f is given by rj b f — S 2 /(S 2 + B 2 ), where S 
is the FWHM source size and B is the FWHM beam size. If all 
transitions included in the analysis are optically thin and rj\,f is 
set to 1 (B « S), the method is very similao to the rotation 
diagram analysis with N mo \ then representing a beam averaged 
column density. 

In selected cases we will also check the results and refine the 
models obtained with the modified rotation diagram analysis by 
employing a more accurate treatment of the line excitation and 
radiative transfer. We have here adopt ed the accelerated lambda 
iterat ion (ALI) technique outlined by Rvbic ki & Hummer! (1 1 99 ll 
1992). The ALI model cloud consists of spherically concentric 
shells and allows only for radial gradients of the physical param- 
eters (kinetic temperature, molecular hydrogen density, molecu- 
lar abundance, and radial velocity field). Moreover, and in con- 
trast to the rotation diagram method, collisional excitation is in- 
cluded in the analysis, so the collision partner (here H2) density 
is a physical input parameter. The code us ed in the present work 
has been tested by Maercker et al.l (120081) and it allows the in- 
clusion of dust as a source of continuum emission in the shells. 
When the radiative transfer has been iteratively solved using the 



The only difference is that in the normal rotation diagram analy- 
sis a straight line is least-square fitted to quantities proportional to the 
logarithm of the line intensities but here the fit is performed directly by 
minimizing the sum of the squared and error-weighted differences of 
observed and modelled line intensities. 
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ALI approach, a model spectrum is produced by convolving the 
velocity dependent intensity distribution of the model cloud with 
a gaussian beam. 

4.1. Analysis of the formaldehyde lines 

As can be seen in Table|4]all four H2CO lines are clearly detected 
in the three sources (D-peak, PI, and CHsOH-peak). The three 
para lines at 218 GHz have all been observed simultaneously so 
their relative strengths are not affected by any pointing or cal- 
ibration uncertainties. The only ortho transition, the 3 12 - 
line at 225 GHz, was observed in a separate frequency setting. In 
our rotation diagram analysis we treat the ortho and para H2CO 
lines together by assuming that the population distribution is de- 
termined by T ml also between the states of different symmetry 
(through some exchange reaction or formation mechanism). For 
H2CO the lowest ortho level ( 1 11 ) is about 15 K above the low- 
est para level (Oo,o)- For D2CO the reverse situation is in effect 
and it is the lowest level Oo,o that is an ortho state while the low- 
est para level 1 y is 8 K higher up in energy. For a very low r rot 
(during molecule formation) most of the molecules will be in 
the lowest energy state (para for H2CO and ortho for D2CO). On 
the other hand, if r rot is much greater than the energy difference 
of the symmetry states, the ortho to para population ratio will 
be governed by the statistical weight ratio which is o/p = 3 for 
H 2 CO and o/p = 2 for D 2 CO. 

In Table [7] the results of the modified rotation diagram anal- 
ysis are shown. Here the number of lines used in the analy- 
sis for each molecule is tabulated together with best fit val- 
ues of r iot , N mo \, and 7/bf. In the last column, we list the min- 
imum and maximum optical depths of the used transitions in 
the analysis. If there are not enough of lines or all lines are 
optically thin, one or two of the parameters have been set to 
a result obtained by another molecule in a previous fit. For in- 
stance, in the D-peak source, all three parameters could be de- 
termined in the analysis of H2CO, while for the optically thin 
D2CO lines, ?7bf had to be set to the value found for H2CO. 
Using the same rjw for all formaldehyde isotopologues will 
also make the determined column densities directly compara- 
ble with each other. Interestingly, the T ml obtained for D2CO, 
17.4 K toward the D-peak is lower than that found for H2CO 
(r^ = 22.5 K). We assume this is an effect of subthermal exci- 
tation and difference in optical depths, where the higher optical 
depths for H2CO make the excitation more efficient via photon 
trapping as compared to D2CO. A lower rotation temperature of 
the less abundant for maldehyde isotopolog ues was also seen in 
IRAS 1629 3 -2422 bvlLoinard et all d200db and in several other 
sources bv lParise et alT i 20061) . For the other two sources (the PI 
and CH3OH peaks) the number of detected D2CO transitions is 
not sufficient to allow for a T rot determination and we assume 
that the ratio of 7 , rot (D2CO)/7' rot (H2CO) is the same in these two 
sources as in the D-peak source. The deduced T mt is higher in the 
PI and CH3OH peaks than in the D-peak. This could be expected 
because the ratio of the 3o,3 - 22,o and 37 7 ~ 2?.i H7CO lines 
is a g ood measure of kinetic temperature dMangum & Woottenl 
119931) and this ratio is highest toward the CH 3 OH-peak. Of 
course, from the models we can also estimate the intensity of 
lines not included in the analysis. In particular, we find for the 
HDCO li,i - Oo,o line that for the determined model parameters 
toward the D-peak its expected intensity is about 5 times below 
the noise level, consistent with our non-detection. 

The derived H2CO rotation temperat ures for the three c ores 
are quite similar to the ones obtained bv lParise et all d2006l) for 
other low-mass protostar sources. Likewise, the derived H2CO 



Table 7. Analysis results for H 2 CO, H' 3 CO, HDCO, D 2 CO, 
CH3OH, SO, and SO2 from using the modified rotation diagram 
technique 



Molecule 


No. of 


T mt 


iVmol 


lit 






lines 


(K) 


-7 

cm 










D-peak (0", -30") 






H 2 CO 


4 


22.5 


2.22 x 10 14 


0.446 


0.11,2.70 


D 2 CO 


6 


17.4 


3.17 x 10 13 


(0.446) 


0.01,0.54 


HDCO 


1 


(17.4) 


2.37 x 10 13 


(0.446) 


0.23 


H' 3 CO 


1 


(17.4) 


4.09 x 10 12 


(0.446) 


0.06 


CH 3 OH 


6 


6.8 


7.21 x 10 14 


1 


0.02,0.62 


SO 


2 


18.3 


2.38 x 10 14 


(1) 


0.50,1.21 


34 SO 


1 


(18.3) 


9.67 x 10 12 


(1) 


0.02 


so 2 


2 


19.9 


4.23 x 10 13 


(1) 


0.01,0.07 






PI (C 


",+30") 






H 2 CO 


4 


25.9 


1.90 x 10 14 


0.469 


0.10,1.91 


D 2 CO 


2 


(20) 


1.66 x 10 12 


(0.469) 


0.01,0.02 


HDCO 


1 


(20) 


6.21 x 10 12 


(0.469) 


0.05 


H' 3 CO 


1 


(20) 


1.71 x 10 12 


(0.469) 


0.02 


CH3OH 


6 


7.4 


1.24 x 10 15 


1 


0.05,1.20 


S0 2 


2 


23.4 


9.91 x 10 13 


(1) 


0.04,0.12 






CH 3 OH-peak(-30",0") 






H 2 CO 


4 


26.3 


2.66 x 10 14 


0.325 


0.14,2.60 


D 2 CO 


2 


(20) 


1.93 x 10 12 


(0.325) 


0.01,0.03 


HDCO 


1 


(20) 


2.38 x 10 12 


(0.325) 


0.02 


H' 3 CO 


1 


(20) 


3.58 x 10 12 


(0.325) 


0.05 


CH3OH 


6 


9.2 


1.26 x 10 15 


1 


0.11,1.46 


so 2 


1 


(20) 


8.33 x 10 12 


(1) 


0.01 






S-peak ( 


-60", +60") 






S0 2 


3 


19.7 


7.63 x 10 14 


1 


0.04,1.21 


34 S0 2 


2 


(19.7) 


5.63 x 10 13 


(1) 


0.04,0.05 



A value surrounded by (...) indicates a fixed parameter. 



column densities are just slightly higher than the corresponding 
rotation diagram values of Pa rise et ail (120061) . The filling factor 
?7bf = 0.446 determined toward the D-peak is in agreement with 
the filling of ~ 0.5 determined from the D2CO 4o,4 - 3o,3 source 
size in Sect. ( I3.lt . 

4.2. Analysis of the methanol lines 

The CH3OH analysis using the modified rotation diagram tech- 
nique is based on the integrated line intensities listed in Table [5] 
In all sources 6 CH3OH lines with lower state energies ranging 
from 23 to 46 K, have been used. In the D-peak source a couple 
of lines are just marginally detected but are included in the best 
fit analysis since they are weighted with the uncertainty and do 
not affect the fit significantly. The line feature at 241904 MHz 
is a blend of two CH3OH lines of about the same energy and 
A-coefficient and we have not been able to clearly resolve them 
into individual components. All lines at 241 GHz were observed 
simultaneously. The only CH3OH line not belonging to the 241 
GHz 5-4 line forest is the 42,2 - 3 12 E-line which is located in 
the 218 GHz H 2 CO-band. 

Only one line from the methanol A-species has been ob- 
served and included in the analysis. Just like in the case of or- 
tho and para formaldehyde, the A- and E-species of methanol 
are treated together. The energy difference between the lowest 
A and E-methanol states is 8 K with the A-species having the 
lowest energy. Furthermore, we only consider methanol to be in 
its lowest torsional state. 

The best fit results in the CH3OH analysis have been entered 
in Table|7] In all three sources we find a low r rot of 7 to 9 K. Also 
the 77bf is found to be very close to 1 . The much lower rotation 
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temperatures found in the CH3OH analysis as compared to the 
H2CO and D2CO results suggest that the excitation of CH3OH 
is quite sub-thermal and a more elaborate treatment of the exci- 
tation and radiative tr ansfer is needed (as was demonstrated by 
Bach iller et al.l dl998l) in their Fig. 9a). This is also supported by 
the fact that none of the CH3OH fits are good since several of the 
modelled line intensities deviate substantially from the observed 
line intensities. 



4.3. ALI modelling 

Our aim with the non-LTE modelling using the ALI code is 
to construct a homogeneous (in terms of physical parameters 
but not in terms of excitation which may vary radially) spher- 
ical model cloud, for each of the source positions (D-, PI, and 
CH 3 OH-peaks), that produces the observed H 2 CO and CH 3 OH 
spectra. The results from the modifed rotation diagram analy- 
sis above are used as a guide when adopting the model cloud 
parameters. The ALI setup for the statistical equilibrium equa- 
tions includes collisional excitation rates and for formald ehyde 
we use the He-H2CO collisional coefficients calculated bv lGreenl 
( 1991). The coefficients have been multiplied by 2.2 to approx- 
imate the H2-H2CO collision syst em. For the H2-CH3 OH colli- 
sion coefficients we adopt those of Pottage et al. (2004). For both 
species only levels below 200 K are included and we divide the 
model cloud into 29 shells. We do not include any radial velocity 
field in the source and instead we simply use a microturbulent 
velocity width, v tur b, that reproduces the observed line widths. 
For an optically thin line, the FWHM line width is V41n2 v^rb- 
The results of the three ALI models have been summarized 
in Table[8] The listed cloud radii R (in the first column) are based 
on the corresponding H2CO beam filling results (Table 13 and a 
distance to p Oph A of 120 pc. The tabulated column densities 
are the source-averaged values 2 N pea ^/3, where /V pea k is the col- 
umn density through the center of the spherical model cloud. It 
should be noted that the derived /^(FL) (and hence A^Fb)) is a 
compromise value since, for all three models, a 20% lower value 
yields better fits for CH3OH while a similarly higher value is 
optimum for H2CO. However, as the differences are less than 
about 50% any real difference in n(H2) is within the uncertain- 
ties of the adopted collisional coefficients. The ortho-to-para ra- 
tio is about 2 for all three models but since the ortho results are 
based on a single optically thick line (r ~ 4) only and not ob- 
served simultaneously with the three para lines, we cannot really 
exclude ortho-to-para ratios in the range 1-3. The derived H2CO 
column densities are about 30% lower than those derived from 
the modified rotation diagram analysis. In contrast to A^FLCO) 
and A^FL) which do not vary much among the three cores, the 
CH3OH column densities vary from 5 x 10 13 cm~ 2 in the D-peak 
to 4 x 10 14 cirT 2 at the CF^OH-peak position. In addition, the 
E/A CH3OH ratios seem to be close to 1 and since the A-line 
is observed simultaneously with the E-lines the determined A/E- 
ratio of 1 is less uncertain than the estimated ortho/para ratios for 
formaldehyde. The optical depth for the strongest of the CH3OH 
lines are about 0.5-0.6 toward the CFbOH-peak. The total H2 
mass of a core is about 0.1 - 0.2 M . 

To our knowledge, there are no collision coefficients avail- 
able for the collision system H2-D2CO, so we cannot use our 
ALI model for D2CO. However, using the H2-H2CO collision 
rates we check what happens to the excitation temperatures 
when adopting a factor of 7 lower column density (the ratio of 
7V(H 2 CO) and 7V(D 2 CO) in Table |7]is close to 7). We then find 
that the excitation temperatures of the 3 1.2 - 2 1,1 transition drop 



between 20-30% at different radii. This drop in excitation tem- 
perature is in line with the 20% lower r, ot = 17.4 K found for 
D2CO as compared to H2CO value of T rot - 22.5 K found to- 
ward the D-peak. Hence, it is quite plausible that the lower ro- 
tation temperature found for D2CO as compared to H2CO can 
be explained by less efficient photon trapping in the excitation 
of D2CO. We also checked the influence on our results by con- 
tinuum emission of dust, and the inclusion of a dust component 
with standard dust parameters (7d = 22 K, gas-to-dust mass ratio 
of 100) had negl i gible impact on our modelling results. 

iLiseau et al.l d2003l) report results from a large velocity gra- 
dient (LVG) analysis of the p Oph A using CH3OH 2 - 1 and 
3-2 line forest data, at 96 and 145 GHz, obtained by SEST at 
an angular resolution larger (52" and 35") than for the present 
study. The analysed data were averaged from spectra spaced by 
30" in the north-south direction around the P2 position (see 
Fig. [TJ and their spectra is thus a partial blend of the contri- 
bution of the three cores. This is especially the case for their 
2-1 observations. However, their LVG results (7\i n = 20 K, 
n(H 2 ) = 4.5 x 10 5 cm" 3 , X[CH 3 OH] = 2.7 x 10~ 9 ) are in good 
agreement with those reported here. 

4.4. Analysis of other molecules 

In the S-peak position we have detected three SO2 lines and a 
modified rotation diagram analysis gave best-fit results for ^bf 
very close to 1 . For SO2 we find a rotation temperature of about 
20 K and a column density of N(S0 2 ) = 7.6 x 10 14 cnT 2 . The 
same analysis was made for two of the three detected 34 S02 lines 
and we get 7V( 34 S02) = 5.6 x 10 13 crrT 2 when using the rotation 
temperature of the more common variant. The lli,n - 10o,io 
34 S02 line was excluded in the fit since its observed strength 
is incompatible (too strong by about a factor of 2) with the ro- 
tation temperature found for SO2, cf. Table [6] The reason for 
this is unclear but could be due to an unknown blend or a non- 
LTE effect involving K a = states. Likewise, the two SO lines 
detected toward the D-peak position result in a rotation temper- 
ature of 18 K and a column density of 2.4 x 10 14 cirT 2 (again 
assuming ?/bf = 1). In the PI position, the two detected SO2 
lines give T mt — 23 K and a beam averaged column density of 
9.9 x 10 13 cm -2 . Furthermore, assuming r, ot = 20 K we arrive at 
a beam averaged SO2 column density of 8.3 x 10 12 crrT 2 toward 
the CH30H-peak. The SO and SO2 results have been entered in 
Table 

The analysis results of the other molecules are summarized 
in Table [9] The column densities listed here are all calculated 
under the assumption of optically thin emis sion. The C 18 Q(3-2) 
column densities are based on the data from lLiseau et al. 
and have been derived using the kinetic temperatures derived 
from the ALI analysis as excitation temperature as the lower 
C 18 transitions are expected to be thermalized at the high densi- 
ties of the cores. For the S-peak, the SO2 rotation temperature of 
20 K was used as excitation temperature for all molecules. Also 
tabulated is the N(H2) column density deduced from 7V(C 18 0) 
and assuming a C ls O abundance of X[C I8 0] = 2 x 10~ 7 (stan- 
dard value of unde pleted gas interst e llar g as, iFrerking et al.l 
(1982), but see also Wouterl oot et al.l d2005l) for a discussion 
applicable to the p Oph cloud). In this context it is enlight- 
ning to estimate how much the cores contribute to the observed 
C 18 0(3-2) integrated line intensity. Adopting the physical pa- 
rameters (Table [8) and assuming X[C 18 0] = 2 x 10~ 7 we find, 
using the ALI code for C 18 0, that the cores make up 31%, 
37%, and 67% of the observed C ls O(3-2) emission in the D- 
, P1-, and CH30H-cores, respectively. Firstly, this tells us that 



9 



P. Bergman et al.: Deuterated formaldehyde in p Ophiuchi A 



Table 8. Model cloud properties and results for the ALI analysis 



Source 


R/I0' b 


v turb 




n(H 2 ) 


MH 2 ) 


M(H 2 ) 


X[H 2 CO]xl0" 


N(\l 2 CO) 


X[CH 3 OH]xl0 y 


JV(CH 3 OH) 




(cm) 


(kms" 1 ) 


(K) 


(cm" 3 ) 


(cm" 2 ) 


(M ) 


para 


ortho 


(cnT 2 ) 


A 


E 


(cm" 2 ) 


D-peak 


3.3 


0.4 


24 


6 x 10 5 


3 x 10 22 


0.16 


1.8 


3.6 


1.5 x 10 14 


0.9 


0.9 


5.0 x 10" 


PI 


3.4 


0.6 


27 


7 x 10 5 


3 x 10 22 


0.18 


1.7 


3.3 


1.5 x 10 14 


2.0 


2.0 


1.2 x 10 14 


CH 3 OH-peak 


2.7 


0.6 


30 


1 x 10 6 


4 x 10 22 


0.14 


1.7 


3.3 


1.8 x 10 14 


5.3 


5.8 


4.0 x 10 14 



the C 18 abundance in the cores cannot be much higher than 
2 x 10~ 7 as they then would produce too much C is O(3-2) emis- 
sion. Secondly, an appreciable part of the observed C I8 0(3-2) 
emission is likely to arise in a lower H2 density 10 5 cm~ 3 ) 
environment. This would explain the higher N(}\2) column den- 
sities deduced from optically thin C 18 (Table |9]l as compared 
to the ALI results (Tab le [8). A similar finding was made by 
IWouterloot et al.l (12005) in their CO study of other regions in the 
p Oph cloud where the cold (~ 10 K) and dense (^ 10 5 cm -3 ) 
cores were surrounded by a warmer (~ 30 K) and less dense (~ 
10 4 cirT 3 ) envelope. Also, given the observed rati o of C 18 Q(3- 



Table 9. Analysis results for other molecules assuming optically 
thin emission 



2) and 13 C 18 0(3-2) of about 23 toward the D-peak dLiseau et alj 
2010), it is likely that the C 18 0(3-2) emission is somewhat op- 
tically thick (t w 2) and hence the listed N(C li O) and N(H 2 ) 
(Table |9]l would be about a factor r/(l - e~ T ) « 2.3 higher. 
The full H2 column density toward the D-peak would then be 
— 1.5 x 10 23 cm -2 . Likewise, the SO column density in 
the D-peak position of 1.1 x 10 14 cirT 2 is lower than the cor- 
responding value of 2.3 x 10 14 cirT 2 (Table [7]) where the latter 
value includes compensation by optical depth. Using the ALI 
model for the D-peak core (Table [8]) we find that adopting an 
abundance of X[SO] = 1.6 x 10~ 8 results in integrated line in- 
tensites for the SO 5e - 4s and 65 - 54 transitions that are close 
to the observed values (Table[6]i. Hence, using the previously de- 
rived physical properties of the D-peak we can also explain the 
observed SO emission. We assume that the bulk of SO emission 
originates in the D-peak core itself and not from the low-density 
envelope. With the same assumption for the optically thin SO2 
lines toward the D-peak, we find an SO2 abundance of X[S02] = 
N[S0 2 ]/N[U 2 ] ~ (4.23 x 10 13 /0.446)/3 x 10 22 ~ 3 x 10~ 9 us- 
ing the beam averaged SO2 column density of 4.23 x 10 13 cm -2 , 
D-peak filling factor of 0.446 (Table [7]) and the D-peak core 
H 2 column density (Table©. Similarly, for the PI and CH3OH 
peak positions we obtain X[S02] = 7 x 10~ 9 and 6 x 10~ 10 , re- 
spectively. Here the value for the PI position is uncertain since 
it is likely that the SO2 emission do not originate in the same 
source as H2CO since they exhibit different emission velocities 
(cf. Table [3]). Using the PI H 2 column density from C 18 0(3-2) 
(Table [9} instead will lower the abundance by about a factor of 
five. Taken together, the SO2 abundances in the D-peak, PI, and 
CH3OH cores seem to fall in the range (0.6 - 3) x 10 9 . However, 
these SO2 abundances are all lower than the S-peak abundance, 
which can be estimated to X[S0 2 ] = 2 x 10~ 8 . 



4.5. Column density ratios 

To estimate the ratios of different formaldehyde isotopologues 
we use the results obtained by the modified rotation diagram 
technique (Table |7). Even for the main isotopic species we 
use the rotation diagram analysis results and not the ALI re- 
sults, since effects like adopting a spherical source will affect 
the results and we rather be consequent in as many aspects as 
possible when estimating column density ratios. The four col- 
umn density ratios; HDCO/H 2 CO, D 2 CO/HDCO, D2CO/H2CO, 
H2CO/HPCO have been calculated for each source position and 



jV(C 18 0) 


N(H 2 y 


iVf 4 SO) 


jV(SO) 


JV(N 2 D + ) 


(cm- 2 ) 


(cm- 2 ) 


(cm" 2 ) 


(cm- 2 ) 


(cm- 2 ) 



D-peak (0", -30") 

1.3 xlO 16 6xl0 22 6.4 xlO 12 1.1 x 10 14 2.3 x 10 12 

PI (0",+30") 

1.4 x 10 16 7 xlO 22 1.6 x 10 13 2.0 xlO 14 2.6x10" 

CH 3 OH-peak (-30", 0") 
7.7 xlO 15 4 xlO 22 <2.2xl0 12 1.2 xlO 14 1.2x10" 
S-peak (-60", +60")* 

7.5 xlO 15 4 xlO 22 5.5 xlO 13 2.6 xlO 14 <4.8xl0 10 



" Assuming JC[C ls O] = 2 x 10"' 
h Adopting T ex = 20 K 



entered into TablefTUl The listed ratio uncertainties have been es- 
timated by including an absolute calibration uncertainty of 10% 
together with the lcr uncertainty due to noise, where the latter 
is the most prominent source of error for the weak lines. The 
HDCO/H2CO ratio of 0.1 07 + 0.015 toward the D-peak is simi- 
lar to the values found bv lParise et al.l d2006l) for low-mass pro- 
testers. The D2CO/HDCO ratio of 1.34 ± 0.19 is higher than the 
Parise et al. (2006) values which all fall in the range 0.3-0.9. In 
fact, it is the only reported D2CO/HDCO ratio that is signifi- 
cantly greater than 1 . 
The quantity 



F = 



(HDCO/H2CO) 2 HDCO/H2CO 



(1) 



D2CO/H2CO D2CO/HDCO 

is also tabulated in T able [TOl This quantity was introduced by 
Rodgers & Charnlevl J2002) to discern between gas phase and 
grain surface formation mechanisms of deuterated formalde- 
hyde. Our D-p eak F-value of 0. 8 is si milar to the values 0.03- 
0.2 found by iRoberts & Millarl d2007l) . In the last column of 
Table [TO] we also list the H2CO/CH3OH column density ratios 
base d on the ALI mode l results. 

iBensch et al.1 d200ll) used observations of rare CO isotopo- 
logues to determine the 12 C/ 13 C ratio toward p Oph C to 65 ± 1 1 . 
Although the p Oph C cloud core is quite far from the p Oph A 
position they belong to the same nearby cloud complex and our 
^CO/H^CO ratios (TablefTUli are in good agreement with their 
value. Consequently, we will adopt this I2 C/ 13 C ratio of 65 also 
for p Oph A and the opacity for the C 18 0(3-2) line toward the D- 
peak can be estimated from th e optical depth of the I3 C 18 0(3-2) 
line of 0.03 dLiseau et alJ2010h as -rfC 18 0(3 - 2)1 * 65x0.03 ~ 
2. This is then consistent with the discussion in the previous sec- 
tion on the C 18 0(3-2) opacity corrections toward the D-peak. 

In Table HQ the SO/ 34 SO and S0 2 / 34 S0 2 column density 
ratios have been listed. As noted in the table, some of the 
ratios are affected by opacity in the main isotopologue and 
shoul d be regarded as low er limits. The interstellar 32 S/ 34 S 
ratio (Lucas & Liszj 19981) appea rs to be close to its solar 
value of 22 dAsplund et al.l 120091) . We can use this isotopo- 
logue ratio to estimate the SO abundances from the expres- 
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Table 10. Formaldehyde column density ratios 



Source 


HDCO/H 2 CO 


D,CO/HDCO 


D 2 CO/H 2 CO 


h 2 co/h; 


CO 


F 


H2CO/CH3OH 


D-peak(0",-30") 


0.107 + 0.015 


1.34 ±0.19 


0.143 ±0.020 


54: 


1 13 


0.080 ± 0.016 


3.0 ±0.4 


PI (0",+30") 


0.033 ± 0.006 


0.27 ± 0.08 


0.0087 ± 0.0024 


111 : 


h32 


0.12 ±0.04 


1.2 ±0.2 


CH 3 OH-peak(-30",0") 


0.009 ± 0.004 


0.8 ±0.5 


0.0073 ± 0.0025 


74: 


t 18 




0.45 ± 0.06 



Table 11. Other column density ratios 



N[H 2 ] gas /N[H 2 ] dust 



Source 


so/ 34 so so 2 / J4 so 2 


N 2 D + /N 2 H + "" 


D-peak 


25 ±5 


0.078 ±0.011 


PI 


13 ±2* 


0.0083 ± 0.0024 


CH 3 OH-peak 


> 18'' 


0.024 ± 0.010 


S-peak 


4.7 ± 0.7* 14 ± 3 





a A'(N 2 H + ) estimated from Andre et al. 
b Not corrected for opacity 
c 3cr limit 



sion 22 Af[ 34 SO]/./V[H 2 ]. For the D-peak we already have an 
abundance estimate, made in the previous section, of X[SO] = 
1.6 x 10" 8 . In the PI and CH 3 OH peaks we then get 2 x 10~ 8 
and 3 x 10~ 9 , respectively. The PI value is quite uncertain, as 
the SO emission velocities differ here from that of H 2 CO (cf. the 
case of SO2 in the previous section), will be lower by a factor 
5 when using beam averaged column densities. Using the val- 
ues in Table [9] for the S-peak we find that the SO abundance is 
X[SO] = 22 x 5.5 x 10 13 /4 x 10 22 ~ 3 x 1(T 8 . 

The ratio N 2 D + /N 2 H + is also tabulated in Table [TT1 and here 
we us e the N2H + (l-0) integrated intensity data from Andre' et all 
(2007) which have a beam size of 26", i.e. almost identical to the 
resolution of ourN 2 D + (3-2) data. The listed jV(N 2 D + )/JV(N 2 H + ) 
ratios have not been compensated for optical depth of the N 2 H + 
lines. 



5. Discussion 

5.1. Interaction with the molecular outflow? 

As described in Sect. ( I3.3b . the observed line profiles for SO and 
S0 2 show broader wings in the north-western corner of our map 
(Pig . [9j» . Also the H 2 CO spectra exhibit line wings in this region 
(Fig. |2]l. The CO red wing emission from the prominent outflow 
emanating from VLA 1623 extends in this direction while the 
blue wing emission is e xtended both in the NW and SE direc- 
tions d Andre et al.| [T990). The SO and SQ 2 emission peaks very 
close to the IR reflection nebula GSS30 dCastelaz et al.l fl985) 
which is at offsets (-79", +28") in our maps. The SO and S0 2 
abundances (3 x 10~ 8 and 2 x 10~ 8 , respectively) are here clearly 
higher than in other positions. The SO and S0 2 line wings could 
very well be associated with GSS30 rather than the VLA 1623 
CO outflow. The CH 3 OH-peak at (-30", 0"), on the other hand, 
could be a result of an interaction with the outflow and a dense 
clump (cf. [j^rgensen et al. 2004). The methanol abundance is 
about 5 times higher here compared to its abundance in the 
D-peak and PI cores. It appears also to be somewhat warmer 
(~ 30 K) than the other cores. 

5.2. Depletion 

The D-peak position where the deuterium fractionation is high- 
est coincides with the emission peak (SMI) in the 1.3 mm 
continuum map of Mot te et al.l (|1998l) . The 1.3 mm flux den- 
sity, S i.3 mm , can be converted into an H 2 column density using 




Aa [arcsec] 

Fig. 11. Gas-to-dust ratio map of p Oph A calculated as the 
H 2 column density from C ls O relative to the H 2 column den- 
sity fr om dust. The gas c olumn data are based on the C 18 0(3-2) 
map of | Liseau et aI1(l2010l) and the dust continuum data are from 
Mott e et al.ldl998lh The gas column data have not been corrected 
for opacity. 



JV(H 2 ) = 4 x lO^cnr^Somm/mJybeanr 1 ) 

assuming 7/dust — 

20 K, a dust mass opacity of K]p, m m = 0.5 cm 2 g" 1 , and a gas- 
to-dust mass ratio of 100 (e.g. Andre et afl [l996) over the entire 



map (see Eq. (1') in iMotte et al. 1998b . This can be compared 
to the H? column density derived from the C 18 0(3-2) map of 
iLiseau et al.l (HblO). Following the discussion in Sect. ( 14.4b we 
here assume X[C I8 0] = 2 x 10~ 7 and an excitation temperature 
of 20 K in all positions. The ratio, N(H 2 ) B!is /jV(H 2 ) dust , of these 
two estimates of A^(H 2 ) are displayed in Fig.Qj] The lowest val- 
ues of this ratio are around 0. 15 and they are found in the central 
part of the cloud. Opacity corrections for C I8 in this central 
part are of the order 2 (which will increase the ratio to 0.3), see 
Sect. (14.41 1 and probably less (closer to 1) further out. This would 
still mean that the A^(H 2 ) gas estimate is a factor of 3 lower than 
that from dust in the central portions of p Oph A. Here the major 
uncertainty is likely to be related to the adopted C 18 abundance 
and dust mass opacity, k\ .3 mm , while influences from temperature 
changes are smaller. The low ratios seen could be due to deple- 
tion of CO and would then correspond to a CO depletion factor 
of 2-3. Given the uncertainty of the adopted dust parameters and 
degree of depletion (ie. C 18 abundance) we estimate that the 
full H 2 column density toward the D-peak, cf. Sect. (14.41 1. is in 
the range (1 - 4) x 10 23 crrT 2 . 

The freeze-out of CO on grains (together with low tem- 
peratures) is thought to increase the deuterium fractionation. 
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Although the D-peak at (0", -30") is located in the part where 
the jV(H 2 ) gas /W(H 2 ) du ' st opacity corrected ratio is low ~ 0.3 - 0.4, 
we do not see any signs that CO depletion is in effect more at 
(0", -30") than in any of the other central positions. Hence, the 
amount of CO depletion (at least as measured by the C ls O emis- 
sion and the 1.3 mm continuum data) is not directly responsible 
to the elevated deuterium fractionation seen in formaldehyde to- 
ward the D-peak. However, the D-peak is located within a larger 
region that is likely to be depleted and this may have been im- 
portant in an earlier stage of the evolution of the D-peak source. 

Interestingly, there is a minimum in the Af(H 2 ) gas /Af(H 2 ) dust 
ratio in the SE part of the map (Fig. ITTb . This minimum co- 
incides with the secondary N 2 D + peak (N6). This would be in 
line with the findings that the amount of CO depl etion and the 
N 2 D + /N 2 H + ratio are correlated in starless cores dCrapsi et alj 



Moreover, the Class source VLA 1623 is also associ- 
ated with a minimum in the A^(H 2 ) gas /A^(H 2 ) dust ratio and possi- 
bly also with a higher N 2 D + /N 2 H + ratio, compare Fig. [10] and 
Fig. Q~T] A correlation b etween CO depletion and N 2 D + /N 2 H + 
ratio has been found by Emprechtinge r et alj d2009l) toward sev- 
eral Class objects. 



5.3. Gas-phase chemistry models 

Pure gas-phase chemistry can lead to elevated deuterium frac- 
tionations for H 2 CO. We have here adopte d the low-metalicity 
(or depleted) model of iRoueff et al.1 (120071) as their warm-core 
(and un-depleted) model shows low deuterium fractionations. 
The steady-state mo del work presented h ere includes the updates 
and modifications of iParise et al.1 d2009l) . and we have used it to 
predict relevant D/H ratios for H 2 CO and N 2 H + . The models are 
shown in Fig.[J_2] The different steady-state D/H model ratios are 
shown as a function of temperature and have been calculated for 
two different H 2 densities; 5 x 10 5 cirT 3 and 1 x 10 6 cnT 3 . The 
modelled ratios show little variation due to the used H 2 density. 
The observed column density ratios (taken from Table [10] and 
Table fTTT) are represented as coloured boxes with size accord- 
ing to estimated uncertainty. The dashed rectangle indicates the 
N 2 D + /N 2 H + ratio toward the D-peak with an opacity correction. 

It is noteworthy that the gas-phase depleted models yield 
F-values, see Eq. dT), very close to 0.5 for T 30 K for 
formaldehyde which are significantly higher than the observed 
values F « 0.1, cf. Table [TOl The absolute H 2 CO abundances, 
(3 - 6) x 10~ n , for these models are much lower than the ob- 
served H 2 CO abundances ~ 5 x 10~ 9 . 



5.4. Grain chemistry models 

Our observations point to the fact that D 2 CO is more abundant 
than HDCO toward the D-peak. This is, to our knowledge, the 
first source where such case is observed. We note that HDCO and 
H?CO showed a similar anomaly towards LI 527 (Parise et al. 
2006), but even in that case, D 2 CO was less abundant than 
HDCO. This anomaly is difficult to explain in terms of pure gas- 
phase chemistry, as discussed above. 

We discuss here the possibility that this anomaly stems from 
grain chemistry. Simple grain models which only account for 
H and D additions to CO to form formaldehyde and methanol 
cannot account for enrichment of D 2 CO compared to the sta- 
tistically expected value. However, including the abstraction 
or exchange reactions of t he type HDCO + H — > DCO + H 2 , 
iRodgers & Charnlevl d2002l) argue that the enrichment is any- 
way limited to HDCO/H 2 CO > D 2 CO/HDCO (or F > 1, see 
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Fig. 12. Gas-ph ase steady-st a te mo dels based on the reac- 
tion network by IRoueff et al.l d2007l) . Shown are the D/H pre- 
dicted ratios; HDCO/H 2 CO (top), D 2 CO/HDCO (middle), and 
N 2 D + /N 2 H + (bottom) as function of temperature and for two dif- 
ferent H 2 densities (5 x 10 5 cirT 3 and 1 x 10 6 cm -3 ). The coloured 
boxes correspond to the observed ratios for the three cores as in- 
dicated in the top panel. The size of the rectangles represents the 
estimated uncertainties. 



Eq. dT])). Our observational resul ts (with F sa 0.1) as we ll as 
those of IParise etaT] d2006h and iRoberts & Millaif d2007h . are 
contradictory to this model prediction, showing that even this 
detailed model may be lacking some processes. 

Some laboratory experiments have aimed in the last years to 
explain the large fractionations of formaldehyde and methanol 
observed towards low-mass protostars . Different types of ex- 
periments have been done (see IWatanabe & Ko uchi 2008], for 
an overview). iNagaoka et alJ (l2005h have exposed CO ices 
with H and D atoms (with D/H of 0.1) and have shown that 
H 2 CO and CH3OH (as well as their deuterated counterparts) 
are formed. In this experiment, they can reproduce the frac- 
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tionation of all deuterated methanol molecu les observed to- 
wards IRAS 16293-2442 dParise et al.l l2004l) . However, they 
never observe an enrichment of D2CO higher than HDCO 
(Naoki Watanabe, priv . com m.). In a second type of ex- 
periment, Hid aka et al.l d2009l) aimed at clarifying the differ- 
ent formation routes by exposing an H2CO ice sample with 
D atoms. Exchanges are observed, showing that deuterium 
fractionation occurs very efficiently along the reaction path 
H 2 CO -> HDCO -> D 2 CO. In this case, the observed D 2 CO 
can become more abundant than HDCO (see their Figure 3). 
Although a direct extrapolation from these experiments is dif- 
ficult, our observations may be the definitive evidence that ab- 
straction and exchange reactions are playing an important role 
in grain chemistry. It is not clear if reproducing our observations 
will requ ire an increased ato mic D/H ratio (from the 0. 1 value 
used by Nagaoka et al. 2005) incoming on the grains, or if the 
longer timescales involved in the ISM chemistry compared to 
the laboratory experiments also can play a role. Detailed mod- 
elling of the complex processes taking place on the grains would 
be needed to definitely settle this question. 

There is no known central source towards the D-peak. So, 
unlike the case of IRAS 16293-2422, there is no protostar that 
can be responsible for releasing the deuterated material into the 
gas. We there fore need to invoke hea ting of the grains by cos- 
mic rays (e.g. lHerbst & Cupp en 2006) or through the release of 
the fo rmation energy of the newly formed species dGarrod et alj 
l2006h . Both these mechanisms can be of importance for desorp- 
tion of grain mantles in starless cores. However, even if these 
mechanisms are responsible of releasing the deuterated material 
into the gas, there is no obvious reason why they would be more 
effective in the SMI core as compared to the other cores. Hence, 
the underlying question why the D-peak shows such a high deu- 
terium fractionation remains unanswered. 

6. Conclusions 

We have observed a very high degree of deuteration of formalde- 
hyde towards a core (SMI) in p Oph A. In this D-peak, 
the deuterium fractionation manifests itself by a very high 
D2CO/HDCO ratio of 1.34 ± 0.19 while the ratio HDCO/H 2 CO 
is 0. 107 ±0.015. In the other parts of this cloud core the degree of 
deuteration in formaldehyde is much lower (although still very 
high compared to the ISM D/H ratio). For instance, at the PI 
position, about V (or 0.035 pc) north of the D-peak, the corre- 
sponding ratios are 0.27 ± 0.08 and 0.033 + 0.006, respectively. 
A similar decrease in deuterium fractionation between the two 
positions is also seen for N2D + /N2H + . The H2CO abundance rel- 
ative to H2 is estimated to be around 5 x 10~ 9 over the central 
core. 

The CH3OH distribution is clearly different from that of 
H2CO and it has its maximum about 45" to the northwest of 
the deuterium peak. The elevated methanol abundance (by about 
a factor of 5 relative to the D-peak) here could be due to an inter- 
action of the outflow with a dense clump. It would be interesting 
to study how the deuterated version s of CH3OH are distributed 
in this source (cf. Parise et al. 20061). 

In order to understand the reason of the very high deutera- 
tion level observed toward the D-peak we have performed gas- 
phase chemistry modelling (for a depleted source). By look- 
ing at the ratio (HDCO/H 2 CO)/(D 2 CO/HDCO) (the F-value, 
see Eq. (Q])) we find that the models result in values around 
0.5 while the observed values are always around 0.1. Also, 
the absolute H2CO abundances obtained in these models are 
too low by a factor of 100. Hence, the gas-phase chemistry 



scheme cannot account for the observed H2CO abundances and 
deuterium ratios. Instead, we advocate that grain chemistry, in 
terms of abstraction and exch ange reactions in th e reaction chain 
H 2 CO -> HDCO -> D 2 CO dHidaka et alJl2009l) . can be respon- 
sible for the very high deuterium fractionations observed in 
p Oph A. However, before being too conclusive about, e.g., the 
required atomic D/H ratio, these grain chemistry model results 
need to be expressed in observable quantities like the H2CO F- 
value. Again, observations of multiply deuterated methanol iso- 
topologues toward the D-peak in p Oph A will provide additional 
insight as they are key ingredients in the grain chemistry scheme. 

Acknowledgements. We are very grateful to Frederique Motte for sending us her 
continuum data of p Oph and to Naoki Watanabe for sending us unpublished lab- 
oratory results. Excellent support from the APEX staff during the observations is 
also greatly appreciated. BP is funded by the Deutsche Forschungsgemeinschaft 
(DFG) under the Emmy Noether project number PA1692/1-1. 



References 

Andre, P., Belloche, A., Motte, F, & Peretto, N. 2007, A&A, 472, 519 
Andre, P., Martin-Pintado, J., Despois, D., & Montmerle, T. 1990, A&A, 236, 
180 

Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 
Andre, P., Ward-Thompson, D., & Motte, F. 1996, A&A, 314, 625 
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Asvany, O., Schlemmer, S., & Gerlich, D. 2004, ApJ, 617, 685 
Bachiller, R., Codella, C, Colomer, F, Liechti, S., & Walmsley, C. M. 1998, 
A&A, 335, 266 

Bacmann, A., Lefloch, B., Ceccarelli, C, et al. 2003, ApJ, 585, L55 
Bensch, F, Pak, I., Wouterloot, J. G. A., Klapper, G., & Winnewisser, G. 2001, 
ApJ, 562, L185 

Castelaz, M. W., Hackwell, J. A., Grasdalen, G. L., Gehrz, R. D., & Gullixson, 

C. 1985, ApJ, 290, 261 
Ceccarelli, C, Castets, A., Loinard, L., Caux, E., & Tielens, A. G. G. M. 1998, 

A&A, 338, L43 

Ceccarelli, C, Loinard, L., Castets, A., et al. 2001, A&A, 372, 998 
Ceccarelli, C, Vastel, C, Tielens, A. G. G. M., et al. 2002, A&A, 381, L17 
Crapsi, A., Caselli, P., Walmsley, C. M., et al. 2005, ApJ, 619, 379 
Di Francesco, J., Andre, P., & Myers, P. C. 2004, ApJ, 617, 425 
Di Francesco, J., Andre, P., & Myers, P. C. 2009, ApJ, 700, 1994 
Emprechtinger, M., Caselli, P., Volgenau, N. H., Stutzki, J., & Wiedner, M. C. 

2009, A&A, 493, 89 
Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590 
Garrod, R., Park, I. H., Caselli, P., & Herbst, E. 2006, Chemical Evolution of the 

Universe, Faraday Discussions, volume 133, 2006, p. 51, 133, 51 
Gerlich, D., Herbst, E., & Roueff, E. 2002, Planet. Space Sci., 50, 1275 
Goldsmith, P. F. & Langer, W. D. 1999, ApJ, 517, 209 
Green, S. 1991, ApJS, 76, 979 

Giisten, R., Nyman, L. A., Schilke, P., et al. 2006, A&A, 454, L13 
Herbst, E., Adams, N. G, Smith, D., & Defrees, D. J. 1987, ApJ, 312, 351 
Herbst, E. & Cuppen, H. M. 2006, Proceedings of the National Academy of 

Science, 103, 12257 
Hidaka, H., Watanabe, M., Kouchi, A., & Watanabe, N. 2009, ApJ, 702, 291 
Jayawardhana, R., Hartmann, L., & Calvet, N. 2001, ApJ, 548, 310 
J0rgensen, J. K., Hogerheijde, M. R., Blake, G. A., et al. 2004, A&A, 415, 1021 
Klein, B., Philipp, S. D., Kramer, I., et al. 2006, A&A, 454, L29 
Larsson, B., Liseau, R., Pagani, L., et al. 2007, A&A, 466, 999 
Linsky, J. L. 2003, Space Science Reviews, 106, 49 
Lis, D. C, Roueff, E., Gerin, M., et al. 2002, ApJ, 571, L55 
Liseau, R., Larsson, B., Bergman, P., et al. 2010, A&A, 510, A98+ 
Liseau, R., Larsson, B., Brandeker, A., et al. 2003, A&A, 402, L73 
Loinard, L., Castets, A., Ceccarelli, C, Caux, E., & Tielens, A. G. G. M. 2001, 

ApJ, 552, L163 

Loinard, L., Castets, A., Ceccarelli, C, et al. 2002, Planet. Space Sci., 50, 1205 
Loinard, L., Castets, A., Ceccarelli, C, et al. 2000, A&A, 359, 1169 
Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodriguez, L. F. 2008, ApJ, 
675, L29 

Lombardi, M., Lada, C. J., & Alves, J. 2008, A&A, 480, 785 
Loren, R. B., Wootten, A., & Wilking, B. A. 1990, ApJ, 365, 269 
Lucas, R. & Liszt, H. 1998, A&A, 337, 246 

Maercker, M., Schoier, F. L., Olofsson, H., Bergman, P., & Ramstedt, S. 2008, 

A&A, 479, 779 
Mangum, J. G. & Wootten, A. 1993, ApJS, 89, 123 
Maret, S., Ceccarelli, C, Caux, E., et al. 2004, A&A, 416, 577 



13 



P. Bergman et al.: Deuterated formaldehyde in p Ophiuchi A 



Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 

Muders, D., Hafok, H., Wyrowski, R, et al. 2006, A&A, 454, L25 

Miiller, H. S. P., Schloder, R, Stutzki, J., & Winnewisser, G. 2005, Journal of 

Molecular Structure, 742, 215 
Miiller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, 

L49 

Nagaoka, A., Watanabe, N., & Kouchi, A. 2005, ApJ, 624, L29 
Nummelin, A., Bergman, P., Hjalmarson, A., et al. 2000, ApJS, 128, 213 
Pardo, J. R., Cernicharo, J., & Serabyn, E. 2001, IEEE Transactions on Antennas 

and Propagation, 49, 1683 
Parise, B., Castets, A., Herbst, E., et al. 2004, A&A, 416, 159 
Parise, B., Ceccarelli, C, Tielens, A. G. G. M., et al. 2006, A&A, 453, 949 
Parise, B., Ceccarelli, C, Tielens, A. G. G. M., et al. 2002, A&A, 393, L49 
Parise, B., Leurini, S., Schilke, P., et al. 2009, A&A, 508, 737 
Pottage, J. T., Flower, D. R., & Davis, S. L. 2004, MNRAS, 352, 39 
Roberts, H. & Millar, T. J. 2007, A&A, 471, 849 
Rodgers, S. D. & Charnley, S. B. 2002, Planet. Space Sci., 50, 1125 
Roueff, E., Lis, D. C, van der Tak, F. F. S., Gerin, M., & Goldsmith, P. F. 2005, 

A&A, 438, 585 
Roueff, E., Parise, B., & Herbst, E. 2007, A&A, 464, 245 
Roueff, E., Tine, S., Coudert, L. H., et al. 2000, A&A, 354, L63 
Rybicki, G. B. & Hummer, D. G. 1991, A&A, 245, 171 
Rybicki, G. B. & Hummer, D. G. 1992, A&A, 262, 209 
Snow, T. P., Destree, J. D., & Welty, D. E. 2008, ApJ, 679, 512 
Turner, B. E. 1990, ApJ, 362, L29 
Turner, B. E. 2001, ApJS, 136, 579 

Vassilev, V., Meledin, D., Lapkin, I., et al. 2008, A&A, 490, 1157 
Ward-Thompson, D., Robson, E. I., Whittet, D. C. B., et al. 1989, MNRAS, 241, 
119 

Watanabe, N. & Kouchi, A. 2008, Progress In Surface Science, 83, 439 
Wouterloot, J. G. A., Brand, J., & Henkel, C. 2005, A&A, 430, 549 



14 



P. Bergman et al.: Deuterated formaldehyde in p Ophiuchi A , Online Material p 1 




Fig. 4. Velocity position diagrams of the H 2 CO(3 ,3 - 2 , 2 ), HDCO(4i, 4 - 3 1>3 ), D 2 CO(4 , 4 - 3 . 3 ), and C 18 Q(3-2) l ines along the 
declination axis at Aa = 0". The molecule is indicated in each panel. The C 18 0(3-2) data are from lLiseau et al 1 <l2010h . The vertical 
scale in each panel indicates the declination offset and the horizontal scale is the vlsr velocity. The intensity scale in all four panels 
is T mh . 
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Fig. 5. Ortho H 2 CO(3i ;2 - 2 W ) and H^ 3 CO(3i, 2 - 2 U ) spectra at declination offsets -30" (bottom, D-peak), 0" (middle) and +30" 
(top, PI). The offsets are relative the center position in Fig.Q] The H^CO spectra (in colour) have been scaled up with a factor of 
10 and appear much more noisy. 
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Fig. 9. SO(56 - 45) and J4 SO(6g - 54) map spectra toward the p Oph A cloud. The J4 SO spectra (in colour) have been multiplied 
with a factor of 5. Scales as in Fig. [2] 



